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ABSTRACT 

We study the two-dimensional kinematics of the Ha-emitting gas in the nearby barred Scd galaxy, 
NGC 6946, in order to determine the pattern speed of the primary m = 2 perturbation mode. The 
pattern speed is a crucial parameter for constraining the internal dynamics, estimating the impact 
velocities of the gravitational perturbation at the resonance radii, and to set up an evolutionary 
scenario for NGC 6946. Our data allows us to derive the best fitting kinematic position angle and 
the geometry of the underlying gaseous disk, which we use to derive the pattern speed using the 
Tremaine- Weinberg method. We find a main pattern speed r2p=22t^ km s^^ kpc^^, but our data 
clearly reveal the presence of an additional pattern speed rip=47l2 kpc~^ in a zone within 

1.25 kpc of the nucleus. Using the epicyclic approximation, we deduce the location of the resonance 
radii and confirm that inside the outer Inner Lindblad Resonance radius of the main oval, a primary 
bar has formed rotating at more than twice the outer pattern speed. We further confirm that a 
nuclear bar has formed inside the Inner Lindblad Resonance radius of the primary bar, coinciding 
with the inner Inner Lindblad Resonance radius of the large-scale m — 2 mode oval. 
Subject headings: Galaxies: spiral - Galaxies: kinematics and dynamics - Galaxies: individual; 
NGC 6946 



1. INTRODUCTION 

Secular evolution of structure in disk galaxies is largely 
governed by internal gravitational processes. In this 
context, density waves play an important role as a driver 
for the internal evolution. Self-gravitating bars contain a 
significant fraction of the disk material and thus have a 
prominent role in the evolution of their host disks. Bars 
are efficient drivers for redistribution of angular momen- 
tum within their host galaxies, leading to mass transfer 
from the outer parts towards the circumnuclear regions 
(e.g., Shlosman, Frank, & Begelman 1989), ultimately 
via structures such as nuclear bars or spirals, continuing 
to the vicinity of their host galaxy's nuclear supermassive 
black hole (e.g., Fathi et al. 2006). An increasing num- 
ber of observational as well as theoretical results suggest 
that secular evolution of disks due to density waves does 
occur, but the details on how the structures and the 
internal dynamics are affected require more inputs, in 
particular to constrain the mass models in the inner 
parts of galaxies where the rotation curve is still rising. 

A key parameter in self-consistent models of barred 
galaxies is the angular rate at which the bar pattern 
rotates, i.e., the pattern speed of the bar, denoted by 
f2p. This parameter is needed to determine whether a 
bar is spontaneous or tidal, to derive the velocity fields 
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as well as mass transfer rates, and hence the densities, in 
perturbed disks, or to establish whether spiral arms and 
any secondary bars are driven by resonance interactions 
with primary bars (e.g., Shlosman & Heller 2002; 
Block et al. 2004). However it has proved complex to 
determine this parameter observationally. 

The most common technique is based on assumed 
knowledge of the locations of the resonances (Buta & 
Combes 1996) which in turn involves understanding the 
behavior of stars and gas at these resonances. Thus, 
purely observationally, ftp can be estimated by the loca- 
tions of the resonances using surface brightness profiles 
of galaxies (e.g., Puerari & Dottori 1997). Alternatively, 
flp can be estimated by matching observations to models 
as suggested by Athanassoula (1992) and successfully 
applied by Salo et al. (1999). A valuable technique for 
the determination of ftp was introduced by Tremaine 
& Weinberg (1984, hereafter referred to as the TW 
method), which allows the measurement of this param- 
eter from kinematic measurements without adopting 
any specific dynamical model. Its application requires 
only measurements of the distribution of intensity and 
velocity of a component that reacts to the density wave. 
It relies on the validity of the continuity equation and 
assumes that the component rotates in a unique and 
well defined pattern. Thus the method has been mostly 
applied to stellar kinematic measurements (e.g., Mer- 
rifield & Kuijken 1995; Corsini, Debattista & Aguerri 
2003; Maciejewski 2006). The TW method is, however, 
in principle applicable for any tracer of a component's 
surface density inferred from the intensity of the tracer 
(Zimmer, Rand, & McGraw 2004). The main uncer- 
tainty of the TW method results is due to the presence 
of dust, although as shown by Gerssen & Debattista 
(2007), these uncertainties do not significantly influence 
the derived pattern speed values. The applicability of 




Fig. 1.— The observed Ha velocity field for NGC 6946. The 
thick diagonal line shows the kinematic major-axis, and the dashed 
ellipses show the location of the corotation radius (outermost), the 
OILR (middle), and the IILR radius of the oval (see section |4]|. 
The dashed lines indicate the arcs inside which the rotation curve 
of Fig. |3] has been averaged, and the parallel dotted lines show the 
position of blocks of 10 apertures along which velocity measure- 
ments have been extracted. The position of the dynamical center 



is set to (x, y) = (0, 0), and the systemic velocity Ksj 



=46 km s 



the TW method to Ha emission was shown convincingly 
by Hernandez et al. (2005) for M 100 and Emsellem 
et al. (2006) for NGC 1068 which also showed some 
evidence of mode coupling between different resonances. 

Here, we apply the TW method to Fabry-Perot obser- 
vations of the Ha-emitting gas in the nearby late-type 
spiral galaxy NGC 6946, with an inclination of 38° 
(Carignan et al. 1990; Zimmer et al. 2004; Boomsma 
2007) and at a distance of 5.5 Mpc (Kennicutt et al. 
2003). NGC 6946 is a grand-design barred spiral galaxy 
which shows three main gravitational distortions: a large 
oval with radius R k, 4.5', a primary bar with ellipticity 
0.15 and radius R « 60", a nuclear bar with ellipticity 
0.4 and radius R ss 8", which is almost perpendicular to 
the primary bar (Elmegreen, Elmegreen, & Montenegro 
1992; Regan & Vogel 1995; Elmegreen, Chromey, & 
Santos 1998; Schinnerer et al. 2006). The disturbed 
spiral-arm morphology has led to the conclusion that an 
m = 3 distortion is superimposed on the primary m = 2 
mode, both with the same corotation radius, R w 4.5' 
(Elmegreen et al. 1992; Crosthwaite 2001). These struc- 
tural properties make NGC 6946 an ideal laboratory 
for a detailed study of the Op, to investigate whether 
the large-scale oval is capable of driving the formation 
and evolution of the inner components by resonance 
interaction. We present the data in section[2l followed by 
the analysis outline and the presentation of the results 
in section [31 and finally the discussion in section H) 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 6946 wa s observed with the FaNTOmM Fabry-Perot 
interferometer (|Gach et al. 1 1200^ mounted on the 1.6 



Fig. 2. — The observed Ho intensity map used in conjunction 
with the velocity field presented in Fig. [T] to calculate the f!p. 
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Fig. 3. — Our derived rotation curve (dots with error bars), the 
average of velocities inside the area outlined by das hed lines in 
Fig. Ul (curves), and the HI rotation curve taken from iBoomsmal 
feOOTi ) marked as crosses. 



meter telescope at the Observatoire du mont Megantic. 
The interference filter, centered at Ac — 6569 A, was 
used to scan the Ha emission-line in 40 channels of 0.21 
A during a total of 120 minutes. The Fabry-Perot was 
tuned to a free spectral range of 391.9 km s~^, over the 
13.7' X 13.7' effective field of view with l'.'6/pix. Complet e 
observation logs can be found in lDaigle et"an (|2006aD . 
The reduction of the data cubes was pe rformed using the 
ADHOCw package (jAmram et al. |[l998f) with consistency 
checks by a pplying the IDL-base d reduction package 
developed bv iDaigle et al. I ()2006bD . Furth er details on 
the re duction procedure can be found in iFathi et al. I 
(|2007l ). After reduction and final cleaning, 29191 
independent velocity measurements remain (see Fig.[T]). 

3. ANALYSIS AND RESULTS 
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In Cartesian coordinates, where the x and y axes are 
aligned parallel to the apparent major and minor axes of 
the galactic disk at inclination, the TW method yields 
the pattern speed 



Q,p sini 



I[x) [V\os{x) - Vsys] dx 
I{x) [x — xq] dx 



(1) 



where x is the position along the major-axis, I{x) is the 
observed intensity, t^os is the observed line-of-sight ve- 
locity and V^ys, the systemic velocity. The nucleus of the 
galaxy is located at xq, and the integrals are performed 
along any cut paral lel to the apparent major-axi s. As 
originally noted by iTremaine &: Weinberg I (jl984f ). the 
use of an erroneous dynamical center could introduce 
errors in the derivation of the f2p, although this could 
be avoided by using an appropri ate weight function to 
elimin ate or reduce the anomalies. iMerrifield fc Kuiiken I 
()1995[ ) refined this technique by normalizing the mea- 
sured quantities with the total intensity in the same aper- 
ture, and obtained 



J_{Vix)l^ 
sini (x) 



(2) 



where {V{x)) is the intensity-weighted 14os-V^ys and {x) 
the intensity-weighted position of the tracer along the 
aperture in question. Thus the calculated flp for each 
aperture does not change, but what changes in equation 
is the position of each aperture in the plot of {V{x)) 
versus (x). As a result, the statistical importance of each 
aperture is evened out, and singularities in positions 
where (x) tends to zero, i.e., near the center of the 
galaxy, are avoided. 

To derive flp, the TW method relies on accurate 
derivation of the geometry of the galactic disk where the 
bar resides. We derive the position of the center, kine- 
matic major-axis, and T^ys by applying the tilted-ring 
method on the observed Ha velocit y field using Tukey 's 
bi-weight mean formahsm ( Mostel lerfc Tukev I Il977f ). 
which is known to be non-sensitive to outliers. Applying 
this formalism has proved to be advantageous for cases 
where outliers, created by bad pixels or strong non- 
circular motions, can bias the derived parameters, and 
thus assigning smaller weights to these pixels delivers 
trustable disk geometry as well as kinematic parameters 
()Fathi et al. II2005D . We derive Vsys^ 46 ± 2 km s~\ 
and the major-axis at 241 ± 2°, both in good agreement 
with the compr ehensive neutral Hydrogen analysis by 
iCarignan et al. I (|l990) and Boomsma (2007). 

For the derived disk geometry and V^ys, the corre- 
sponding rotation curve is illustrated in Fig. [31 where 
cuts along the redshifted and blueshifted sections of 
the disk as well as average rotation curve within the 
bounded section shown in Fig. [1] are overplotted. 
Figure [3] demonstrates that the derived disk geometry 
results in a symmetric cut in the velocity field. Assuming 
a symmetric barred structure and a homogeneous dust 
distribution, a correct dynamical center should imply 
that the rotation curve extracted along the major axis 
is identical to that derived from the tilted-ring method. 
Moreover, the non-circular velocities at the blueshifted 
and red-shifted parts of the disk, should fall and rise with 
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Fig. 4. — The TW method applied to the velocity field presented 
in Fig. [1] by plotting the averaged intensity weighted velocities, 
{V{x)), versus the averaged intensity weighted positions, (x), for 
each three-pixels wide aperture, i.e., 1/lOth of the space between 
the parallel dotted lines in Fig. [T] Open symbols show data points 
from apertures outside the OILR of the large oval, and filled sym- 
bols show the data points inside this radius. 
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Fig. 5. — The angular frequency curve (thick curve) derived 
from the combined Ha rotation curve (inner 9 kpc) and HI rota- 
tion curve (outer kpc) presented in Fig. \3\ Applying the epicyclic 
approximation, we also plot the epicyclic frequency curves, Q — K/m 
for m = 2, 3, 4, and Q + k/2. The pattern speed of the oval, with 
errors, is marked as the horizontal line, and the location of each 
resonance radius is marked with a dot. The two innermost dots 
mark the IILR and OILR radii discussed in section l4l 

respect to the rotation curve in a symmetric fashion. 
This is exactly what we observe in Fig. [31 Furthermore, 
we note that our rotation curv e is in agreement with 
that derived from HI s tudies by ICarignan et al~l ()1990D 
and Boomsma] (|2007[ ). with the advantage that our 
data allows us to resolve the inner regions. 

To derive Op, we set the origin of the coordinate sys- 
tem at the derived dynamical center (marked in Fig. [T]) 
and extract the T^os along successive simulated slits 
parallel to the major-axis, in three pixel intervals. Thus 
the derived velocities within each simulated slit (1/lOth 
of the space between the parallel dotted lines in Fig.[T]) is 
collapsed into one single {V{x)) value, with (x) being its 
intensity- weighted distance along the slit. Alternatively, 
the spectra within each simulated slit could be collapsed 
into one single spectrum from which one single {V{x)) 
value can be derived. Both methods yield similar results. 
Plotting these numbers using equation ([2]), in Fig. [H it 
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TABLE 1 

Different pat tern sp e ed an d c orotation radius values for 
N GC 6946. 1-IPe ton ' (IMT); 2-'Bonnarel ft al. I fTolD; 3- 

IElmegreen e t al. (1998); 4- Crosthwaite (2001); 5- 
izlmmer et al. i (|2004^ ; 6- ischinnerer et al. . (,200611 ; 7- this 

STUDY. 



Method 


Hp (km s ^ kpc 


■^) Corotation (") 


1- Morphology 


60 


78 


2- Morphology 


56 


120 


3- Morphology 




155 


4- Morphology 


66 ±7 


300 ± 60 


5- TW on CO 


39 ±13 




6- Dynamical Model 




270 ± 30 


7- TW, Primary QP 


22+* 


304^1^ 


7- TW, Secondary 


47+4 





is clear that a single pattern speed cannot be assigned 
to all data points. Although thi s trend can also be seen 
in a similar analysis by Zim mer et al. I (j2004l ). these 
authors do not mention or quantify it. This is probably 
due to the small number of points in their {V{x)) versus 
(x) plot (40 data points compared with our 221). We 
investigate this behavior in an iterative way. 

Fitting all the data points, we find an average pat- 
tern speed of 21 km s~^ kpc~^. We use this initial 
value to analyze the angular frequency curve shown in 
Fig. [5l where we find that the outer Inner Lindblad 
Resonance radius (OILR) of the large oval is located 
at around « 50". Using this diagnosis as verified by 
the Kolmogorov-Smirnov test, we find that the points 
for apertures inside this radius (marked as gray points) 
fall on a unique slope clearly different from the outer 
apertures. We thus derive two distinct pattern speeds, 
the primary patter speed, rip=22^^ km s^^ kpc^^, and 
the secondary pattern speed flp—A7+2 km s^^ kpc^^ for 
the region inside the OILR. The errors include errors for 
the derived 14ys and position angle. Using this new value 
for the pattern speed of the oval (compared with 21 
km s^ kpc^ ) does not significantly change the position 
of the OILR, and thus we confirm that the inner 49", i.e., 
the region inside the OILR of the oval rotates at a rate 
different from that of the large oval structure. Our nu- 
merous tests, changing the width of the apertures do not 
induce a significant change in the derived pattern speeds, 
however, we note that for the inner apertures, {V{x)) 
has some contribution from the velocities outside the 
OILR of the oval. Removing all pixels outside the OILR, 
the secondary pattern speed could increase by 50%. 

4. DISCUSSION 

In Fig. Owe find a nu mber of resonances, all confirm - 
mg previous findings (|Elmegreen et al. I Il992l I1998D . 
A direct comparison between our corotation radius 
and previous studies is presented in table [1] showing 
a good agreement bet ween our result and th ose by 
ICrc)sthwaite1 (|2001[ ) and lSchinnerer et al. I (120061) Most 
morphological studies yield underestimates (iPeton I 
Il982t iBonnarel et al. |[T988HElmegreen et al. \[l99^ . 

The disturbed morphology of the spiral arms in 
NGC 6946, has for more than a decade raised the issue 
of the presence of m = 3 and 4 perturbations, all fol- 
lowing the same pattern speed as the main m = 2 mode 



TABLE 2 

The POSITION OF resonance radii WITH ONE-cr ERRORS. 



Corotation radius 304_gg arcsec « 8.30 kpc 

Q - f 165^38 arcsec 4.40 kpc 

O - I 105^^4 arcsec ^ 2.80 kpc 

OILR 49I4 arcsec ^ 1.25 kpc 

IILR 14^^ arcsec ^ 0.38 kpc 



(see ISchinnerer et al. I l2006l . for a detailed discussion). 
Our data confirms the presence of two ultraharmonic 
resonance radii attributed to these higher modes. The 
m = 2 mode, however, induces an inner Inner Lindblad 
Resonance radius (IILR), an OILR, and an Outer Lind- 
blad Resonance radius. The derived OILR radius is of 
the same order of the primary bar, and the IILR radius 
falls just outside the nuclear bar, 7? ss 8" (see table [2]). 

Morphological under-estimation of the corotation 
radius yield an over-estimated value for f2p, by almost 
a factor two. The study tha t best matches our derived 
flp is that bv lZimmer et al. I (|2004D . who apply the TW 
method to CO observations. We confirm that down- 
grading the spatial resolution of our data and deriving 
one unique pattern speed, we get a valu e comparable 
with that derived bv I Zimmer et al. I (|2004f ). Our higher 
resolution data allow a more detailed diagnosis of the 
dynamics of the large-scale oval in NGC 6946, by 
allowing us to derive two distinct f2p values. 

This observational result confirms a "bars within bars 
within bars" scenario, where an evolved large primary 
m = 2 perturbation forms an inner bar whose pattern 
speed could be lar ger b y a fa c tor of up to th ree (e.g^ 
Macieiewski & Soa rke I 120001 : lEnglmaier fc Shlosmanl 
120041 )." The bar' may be found at any position angle 
with respect to the large scale oval, and such systems 
are k nown to form nuclear bars (e.g., lErwin fc Sparkel 
Il999| ). suggesting that multiple bars probably generally 
have different pattern speeds. In Fig. [5l we see that 
flp is consistent with the nuclear bar being located at 
the ILR of t he primary bar in NGC 6 946. Furthermore, 
as found by lElmegreen et al.~l ()1998l ). the nuclear bar 
is almost perpendicular to the larger bar, which they 
interpret as the presence of the well known xi and X2 
orbits of the same structure. Placing the corotation of 
the larger bar at the OILR of the oval agrees with the 
nuclear bar being inside the ILR of the primary bar, 
i.e., the IILR of the oval structure. Such a scenario is 
also known to drive efficient mass flow from the outer 
parts toward the nuclear regions in galaxies (Shlosman, 
Frank , fc Begelman 1989) where a starburs t can occur 
(e.g., iZurita et al.1l2004l; iFathi et al. II2005I) . 

The tight agreement between the sizes of the different 
components with the predicted location of the reso- 
nances confirm that the inner components have formed 
by means of resonance interaction of the main m = 2 
gravitational perturbation, i.e., the large oval and two 
prominent spiral arms. As discussed for example by 
Maciejewski (2004), linear theory is applicable only if the 
gravitational perturbation is not strong. In NGC 6946, 
this is confirmed since the main perturbation is an 
oval distorsion and not a strong bar. The infiow is 
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confirmed by our residual velocity field after subtracting 
the rotational model from the observed velocity field. 
A detailed discussion is beyond the scope of this letter, 
however, the velocity field presented in Fig. [T] shows 
the expected signature of the OILR as a "kink" in the 
contour tracing the systemic velocity. 
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